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Abstract 

The relativistic harmonic oscillator (RHO) model of hadrons is used to study quark stars. The mass- 
radius relationship is obtained and compared with bag model of quark star, using Tolman-Oppenheimer- 
Volkoff equation. In this model, the outward degenerate pressure due to discrete Landau levels and 
Landau degeneracy balances the inward gravitational pressure. Where as in bag model the degenerate 
pressure is due to the standard continuum levels which balances the combined inward pressure due to 
gravitation and bag pressure. So in RHO model, the confinement effect is included in the degenerate 
pressure. We found a qualitative similarity, but quantitative differences in mass-radius relationship of 
quark stars in these two models. Masses and radii are relatively larger and the central energy densities, 
required for stable quark stars, are lower in RHO model than that of bag model. 
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1 Introduction 



The relativistic harmonic oscillator (RHO) model of hadrons was extensively studied and successful to explain 
magnetic moments of baryons [T] , hadron spectroscopy and decays [5] . An attempt was also made to derive 
it from the first principle in Ref. [3]. It is also used to explain nucleon-nucleon force and study nuclear 
properties [Jj. As a further extension of the model, here we apply it to study quark stars. Quark star may 
be a big hadron made up of many quarks, subjected to confinement. Such a matter may be found deep 
inside a neutron star or may be a next fate of a neutron star whenever it's mass exceeds the limit such 
that the degenerate pressure due to neutrons fail to support the gravitational collapse [S]. Such a compact 
stars were studied using different confinement models, different equation of state (EoS) and also models 
with mixture of hadrons and quarks, called hybrid stars [5J. Other exotic quark stars, motivation of these 
studies, and possible candidates for quark stars are reviewed in Ref. [7]. It is a hope that the study of quark 
matter in astrophysics may shed light into perturbative and nonperturbative properties of QCD (quantum 
chromodynamics). It is speculated that certain kind of pulsars may be a possible candidates for quark stars 
[7J. Our aim here is to study nonrotating, cold, pure quark star confined by RHO potential using Tolman- 
Oppenheimer-Volkoff (TOV) equation. Of course, such a study was carried out earlier in Ref. [5] where 
the EoS, needed in TOV equation, was derived by assuming that all quarks are confined to the center of a 
spherical bag by RHO potential. However, we know, from the study of statistical mechanics of particles in 
harmonic oscillator potential [S], that pressure of such a system goes to zero at infinite volume limit. In our 
model, we derive EoS by assuming that quarks are confined at different centers in a finite spherical volume. 
It is similar to the well known problem of diamagntism [3] where electrons are confined at different centers 
in a plane perpendicular to a uniform magnetic field in a finite volume. This model gives finite pressure at 
infinite volume limit, as shown in Ref. [HI EI]- In Ref- we also studied an approximate mass-radius 

relationship of quark stars using Newtonian gravitational theory. 

2 RHO model for quark stars: 

It is speculated by many authors [6j [Jj that, before the massive star collapses to black hole when nuclear 
fuel burns out, it undergoes various stages like white-dwarf, neutron star, hybrid star, quark star or other 
exotic quark stars. Here we are interested in RHO confinement on bulk quark system like quark stars. It 
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is a simplified model of quark star without rotation and hybrid nature of actual quark star. It is a system 
of quarks with overall color neutral, but confined by QCD vacuum, like big hadron. There are such models 
with confinement borrowed from bag models of hadron spectroscopy [10) and studied using both TOV and 
Newtonian formulations. In bag model of quark star, quarks are free inside the bag, but confined by both 
bag pressure and gravitational pressure. In our model, quarks are inside a spherical bag, but each quark 
subjected to a RHO potential with confinement centers anywhere within the bag. This is just the extension 
of the 2-dimensional confinement of electrons in the problem of diamagnetism to 3-dimension. 
From the hadron spectroscopy studies, the single particle energy eigenvalue in RHO is pQ, 



e n = yj(2n + ijc* + ml , (1) 

where n — 1,2, ...oo, C q is the RHO coupling strength and m q , the quark mass. Therefore, quarks inside 
a volume V have discrete Landau levels instead of continuum levels as in the case of free quark (C q = 0) 
inside a bag at infinite volume limit. However, according to Landau's prescription [S], the continuum energy 
levels for C q = have become discrete, but degenerate to accommodate all levels of continuum in RHO. 
This phase-space degeneracy may be written as, 



fin 



(2n + 3) 3 / 2 - (2n + if 2 , (2) 



6tt 2 

a straight forward extension of degeneracy in the problem of diamagnetism [9] . Now following the standard 

procedure of evaluation of degenerate pressure, we first evaluate the energy, 

«f VC 3 np r 

£o = .9/ £ g n e n = P n + 3 ) 3/2 ~ ( 2n + l ) m \ v/(2n+l)C2+m2 , (3) 

n— 1 n— 1 

and the number of quarks, 

gjVCl »* 
6tt 2 

n— 1 n— 1 

where gf is the internal degrees of freedom of quarks and uf is the Fermi level up to which all Landau levels 
are completely filled. Since each level is degenerate, to take account of partially filled level, expression for 
energy may be modified as, 

E » = EI ^T L E [( 2n + 3 ) 3/2 - ( 2n + i) 372 ] y^n +1)^+^2 + (N — N )y/(2n F + 3)C* + m? q , (5) 

n— 1 
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where N is the total number of quarks including the partially filled level. The degenerate pressure is 
P = -■|fa. This degenerate pressure balances the gravitational inward pressure according to the TOV 
equation, 



dP _ M 
df f 2 

and 



[s + P] 



1 



4TTf 3 P 

M 



, { 2M 

f 



-l 



(6) 



— = We , (7) 
where e = energy density, and M is the mass of quark star. All variables above are normalized as, 

As usual, TOV equations are solved, numerically, starting from initial value up, which corresponds to some 
central density, and the radius for which pressure goes to zero is obtained. This is repeated for different central 
densities to obtain the mass-radius relationship, the maximum mass, M max , of the star and corresponding 
radius R are obtained. The actual or unnormalized mass and radius are obtained from 

6^\ 1/2 1 



6^ 1/2 



and plotted in Fig. 1. 

3 Bag model for quark star: 

For comparison, let us briefly review quark stars in bag model [10] . Here the combined pressure due to 
gravitation and bag pressure balances the degenerate pressure to form quark star. Using the standard 
procedure again, the energy, 



p=0 

and the number of quarks, 



Eo = 9f ]T ^p 2 +m\ + BV =%^J o d PP 2 \]p 2 + ™ 2 q + BV 



(11) 



V */£ (12) 
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where pf is the Fermi momentum. From the energy, Eq. degenerate pressure may be obtained and 

using the same TOV equations, Eq. ^ [7|), mass-radius relationship of quark star may be obtained. Now 
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Figure 1: Mass (in terms of solar mass Mq) vs. radius (in km) plot for quark star using RHO model and 
bag model for C q = B 1/4 = 145 MeV and m q = 145 MeV. 



the variables are normalized as 
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Again TOV is solved numerically for different initial values of central densities to obtain mass-radius rela- 
tionship, M and f corresponding to the maximum mass is obtained. From M and f, the actual mass and 
radius are obtained from, 



M = M 



R = R 
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(14) 



(15) 
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which has a similar structure as that of quark star in RHO and mass-radius relationship is plotted in Fig. 
1. C q of RHO is replaced here by B 1 / 4 . Similar mass and radius dependence on bag constant is also 
demonstrated in Ref. [10] using a simple energy principle. 

4 Results: 

The TOV equations, Eq. (0 [7J|, is numerically solved for both the models with different central densities. 
We have considered quark star with two flavors and hence gf — 12. The mass of the star increases with the 



Table 1: Mass, Radius, energy density and quark density of quark star having maximum mass for different 
quark masses in RHO model with C q = 145 MeV. 



m q 


M max 


R 


e c 


n c 


MeV 


M Q 


km 


GeV/fm 3 


fm~ 3 





1.96 


12.55 


0.87 


2.42 


145 


1.7 


10.9 


1.34 


3.22 


nip/3 


1.14 


8.15 


2.14 


4.085 



Table 2: Mass, Radius, energy density and quark density of quark star having maximum mass for different 
quark masses in bag model with B 1 / 4 = 145 MeV. 



m q 


M max 


R 


Sc 


n c 


MeV 


Mq 


km 


GeV/fm 3 


fm~ 3 





1.98 


11.22 


1.05 


2.775 


145 


1.6 


9.4 


1.41 


3.31 


m p /3 


1.02 


6.55 


2.98 


5.32 



central density and reaches a maximum value and decreases as shown in Fig. 1. There is a clear difference 
between the two models. The maximum mass and corresponding radius are tabulated in Tables 1 and 2. 
We have taken C q = B 1 ' 4 = 145 MeV, which is in the range of values used in hadron spectroscopy, and 
obtained the maximum mass for different quark masses, 0, 145 MeV and m p /3, where m p is the proton mass. 
M max in RHO is almost same as that in bag model for massless quarks, but differ for finite masses. M max 
and corresponding radius both are higher in RHO than bag model. M max in bag model is almost same as 
that of Ref. [10]. For massless quarks, quark star mass as well as radius scales as C~ 2 or (B 1 ^ 4 )^ 2 , the 
confinement strengths. Therefore, larger the confinement strength smaller the mass and the radius. From 
the mass-radius plot, Fig. 1, we see that the stable quark stars in RHO model always have larger mass and 
radius compared to the bag model. From Tables 1 and 2, the central energy densities required for stable 
stars are lower in RHO model. For stable stars, mass increases with central energy density, and hence have 
lower central energy density than that of maximum mass star. 

5 Conclusions: 

Here we have applied the RHO model of hadrons to study quark star, using TOV equations and compared 
with extensively studied bag model. In RHO model, confinement leads to discrete Landau levels and as- 
sociated Landau degeneracy, like in the problem of Landau's diamagnetism, instead of continuum levels as 
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in the case of bag model. Hence, the EoS is different from that of bag model and mass-radius relationship 
is obtained for quark star using this new EoS and TOV. On comparison with that of bag model, there is 
a clear difference in mass-radius relationship. For finite quark masses the maximum mass of quark star 
and corresponding radius are higher in RHO than that of bag model because of lower average confinement 
pressure inwards. Masses and radii of stable quark stars are higher, but central energy densities are lower 
in RHO model than that of bag model. 
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